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On Systems of Linear Differential Equations of the 

First Order.* 

By Maxime Bocher. 



It is my object in this paper to consider the system of n differential equa- 
tions : 

^ ='J>y% + A, (t = l, 2, ....n), (1) 

where the coefficients a and fi are functions, not necessarily analytic, of the real 
variable x which, throughout the interval 

(/) a<x<b, 

satisfy the following condition : 

A. The functions a tj and fi t shall have at most a finite number of discontinui- 
ties in (I) , and the integrals 

I | a,;,- \dx, I (3 t dx 

shall converge, when extended over any portion of (I) . 
We add the definition : 

B. A set of n functions of x (y 1} y % , .... y n ) are said to form a solution o/(l) 
when, and only when, they are all of them continuous at every point of (I), and, 
except at most at a finite number of points of this interval, have derivatives which 
satisfy (1). 

The first theorem we wish to prove is the following : 

I. Any set of n constants (cj, <%, . . • • c n ), and any point c of (I) being given, 

there exists one, and only one, solution (y ly y % , .... y n ) of (l) which satisfies the 

initial conditions : 

yi(c)= Ci , (* = 1, 2, .... n). (2) 

* This paper was read before the American Mathematical Society, February 22, 1902. 
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To prove that at least one such solution exists, we use the method of 
successive approximations.* For this purpose consider the series 

y? + yf ] + yf + •••• (»=i, 2, .... n), (3) 

where & = « + fadx, 

and, when &>0, ,_. (." = 1, 2, . . . . n), (4) 

»P , =jT(S«««?" 1) ) <fa - 



i=i 



It is clear that the integrals here written all converge, and represent func- 
tions which are continuous throughout (/) . 

In order to prove the convergence of the series (3), let M be a positive 

constant such that 

\yT\SM, (»=1, 2 n), (5) 

and let <p (x) be a function which has at most a finite number of discontinuities 
in (/) , for which the integral 

/ 4> (x) dx 

converges when extended over any portion of (/), and such that 

KI<<M»0>t (»./=!, 2 n). (6) 

We can now establish the fundamental inequalities : 

n lc \ I <p(x)dx , ,. , „ . 

These inequalities reduce to (5) when & = 0. We shall therefore have estab- 
lished their truth by the method of mathematical induction if, assuming (7) to 

* This method was used for this purpose for the first time by Peano ; cf. Mathematische Annalen, 
vol. 32 (1888), p. 450. Peano, however, restricts himself not merely to homogeneous equations (i. e. 
assumes £i = P 2 = . . . = ]3» =0), but also restricts the coefficients a v to be continuous throughout (I). 

t That such a function exists will be seen from the example : 

i — n j = n 

2 2 K I' 

i = l 4 = 1 
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hold when &<&i, we can prove that it holds when h = \ + 1. Now we have : 

By integrating by parts we find 

4»(a;)<teJ 

j c * ^ v. 9 w "v *" ~ — fc + i ~ • 



(9) 



The substitution of (9) in (8) gives us (7) for the case h = \ + 1 . 
If, now, we write 

N=£<p{x)dx, (10) 

we deduce at once from (7) the further inequality : 

If, then, we compare (3) with the convergent series of positive constant terms, 

"fum, („) 



* = 



we see that the series (3) are absolutely and uniformly convergent throughout 
(I). Writing 

yi = $ ) +y? + yP+ ...., (»= 1, 2, . ... n) (13) 

we see that the functions {yi, y%, • • • • y„) thus defined are continuous through- 
out (/). Moreover, since, when x=c, all the terms except the first of these 
series vanish, while yf ] (c) =c it it follows that y t (c) = c t . 

In order to prove that (y x , .... y n ) form a solution of (1) according to the 
definition (B) , it will then be sufficient if we can prove that, x being any point 
where the coefficients a tj and /3 4 are continuous, (y ly .... y n ) have derivatives 
at % which satisfy (1). 

For this purpose let us surround x by a short interval (J) throughout 
which the coefficients a it and /3< are continuous. We have : 

2 a H y s + fi t = ft + a^ 2/f + a i3 2 3^ + ••••+<£ ^ 

j = l * = fc=0 * = 

= A +2 a ^ 0> +2%^' + • • • • (»' = 1, 2, . . . . n). 

3=1 j=i 
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This last series clearly converges uniformly throughout (J), since the series 
(13) do so. Itis, however, precisely the series which we should obtain by differenti- 
ating (13) term by term. Accordingly y ( has a derivative throughout (J), and this 
derivative is represented by the series last written. Moreover, the formula just 
obtained shows that this derivative satisfies (1). Thus our proof is complete. 

It remains merely to prove that two solutions (y ly .... y n ) and (y 1} .... y n ) 
of (1) which satisfy the conditions 

yi{ G ) = y~i{. G ) (*= 1, 2, n) 

must be identically equal. For this purpose we notice that if 

Vi — yi — 2/7, (* = li 2, n) , 

then (>7 X , *7 3 , .... >7 m ) is a solution of the system of homogeneous linear equations 



dm _ 
dx 



3 — n 



^a-ijVj' 



(i=l, 2, ....n). 



(14) 



^=i 



Our theorem will, therefore, be established if we can prove the following 
special case of it : 

If iyix, viz-, • • • • J7„) is a solution of (14), and if 

Vi(c) = 0, (»=1, 2, .... n), 
then vji (a) = , (i = 1, 2, n) . 

If possible, let h be a point in (/), where all the functions m do not vanish. 
Then we can, by the portion of theorem I already proved, find n — 1 solutions 
of (14), 

Orf°, nt\ •■•;*%) (»■=!, 2, ....n— 1) 

such that the determinant 



£>(*) = 
does not vanish when x = A .* 



>7i 


^ • 


• • M n - 1 ' 


*7 2 


f > • 


• • v^ 


*7« 


^ • • 


• ■ e- J) 



* If, for instance %,(ft) rf= 0, we may determine these solutions by the conditions : 

vP (h) -vP(h) = = v%- P (h) = 1 , 

j#> (ft)=r0 when id^j. 
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Let us now form the derivative of D (x) , writing this derivative as the sum 
of n determinants, each of which is obtained from D by differentiating all the 
elements of one row. By a simple reduction, we find 

D' (x) = (a n + a^ 4- + «»») D (x) . 

We have, therefore, 

D(x) = Ke' / ° , 

where h is a constant.* Accordingly, D(c) = 7c. But we obviously have 
Z>(c) = 0. Accordingly, D(x)= 0. This, however, contradicts the fact above 
mentioned that D(h) qjb 0. The assumption that there is a point where all the 
functions rii do not vanish has thus led us to a contradiction, and thus our theo- 
rem is proved. f 

The second fundamental theorem which we wish to prove states, roughly 
speaking, that a small variation in the coefficients of the system (1) and in the 
initial values c t produces only a small variation in the solution. The theorem 
is this : 

II. The positive constants B, G and e being given, and also a positive func- 
tion $> (x) which has at most a finite number of discontinuities in (I) and for which 

I 4> (x) dx converges when extended over any portion of (I) , then a positive h exists 
such that if 

]Z n (* = 1, 2, . . . . n) 

are any two systems- of equations whose coefficients satisfy condition (A) and also 
(except possibly at a finite number oj points) the inequalities : 

Wi S — ay\<8, \a {J \<<p(x), |ot. t | <^»(a?), 

(*= 1, 2, .... n) 



|A-AI<5, \fPifa\<B, \fj.dx 



<B, 



* Of., for example, Sauvage, " Theorie generale des systemes d 'equations differentielles lineaires et 
homogenes " (1895), p. 13. This thesis was reprinted from the Annales de la Faculte des Sciences de 
Toulouse, vols. 8 and 9. 

t The proof here given is an immediate generalization of the proof of a similar theorem given by 
Sturm (Liouville's Journal, vol. I (1836), p. 109), for a single linear differential equation of the second 
order. 
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(where the last two ineqalities are supposed to hold no matter what two points of (I) 
are taken as t7ie limits of integration), and if c is any point of (I) and (y lf .... y n ) 
and (y x , .... y n ) are solutions of (1) and (T) respectively, for which 

Jy t (c) - y,(o)|<a, \y t (o)\<0, |y,(«)|<0, (.'= 1, 2, . . . . n), 
then i — / \ { m ^ /» = 1, 2, n\ 

M*) — #(*)!<*. U<»<& /* 

To prove this theorem, we express the solutions yi and y { by the method of 
successive approximations, as explained above, in the form : 

(3) y 1 = y? ) + yi 1) + y? ) + ..... 
(3) y i = y? ) + yT ) + yf+ •••-, 



„ _ _ (*== 1, 2, . . . • n). 



Writing 5 -f- C = Jf, we see that every term of (3) and also of (3) is in 
absolute value less than or at most equal to the corresponding term of (12). 

Let us choose a positive integer I, so that the remainder of (12), after the I th 
term, is less than e/3. Then we have : 

\v<-*i\ <5f l^-^l + -T- ' (<= i, a, • • • • »). 

Our theorem will therefore be proved if we can show that S (which is as 
yet wholly unrestricted) can be so chosen that 



i -yi l< 3 r U=0, 1, .... -J- 



\y\ 

Let us assume that, a positive yj being given, a positive 5 exists such that 

|3)_ w W|/„ A* = 1,2, w \ , * 

Then, since we know that such a 8 exists, when k= 0,* we shall have estab- 
lished our theorem by mathematical induction if we can show that by still 
further decreasing $ (if necessary) (15) can be made to hold when h=-m+l. 
For this purpose, let us take S so that 

Yt . m ! 



h< 



2nM{nN) m (b — a) 



* We have then merely to take 6 < y/2, and <J< v/% (6 — a). 
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and then let us still further decrease 8 (if necessary), so that : 

^)_^) [< _^_ > (»=l,2,....n). 



Now we have : 

j = n 



yjP+1) _ j^+D = ^ 2 fog") - «,-#>) tfc 

i=i (t= 1, 2, w). 

= r s p«5^- ^ m> ) + *f° k - °*H ** 

Accordingly, 

|^-+M _ «$«+!> |< I f(nA (x) . -* + nM&^ . ^ZjJg' ) dx\<*i, 

'** ^* '= |«/« V rv ' 2nN ml 2nM(nN) m (b — a)/ I 

(i = 1, 2, . . . . n). 

Thus our theorem is proved. 

As an immediate application of the theorem just proved, we mention the 
following theorem : 

If corresponding to the values k = 1, 2, 3, . . . . , we have an infinite number of 
systems of equations 



T" =2°W + ^ (*•/= I- 2, ... . «), (16) 

i/ iAe coefficients of each of these systems of equations satisfy conditions (A),* awcZ */*, 

in a region consisting of all except a finite number of the points of (I) , affi (x) and 

/?f > (as) approach uniformly the functions <x y and /? s as & becomes infinite, and if the 

constants cf } are such that 

lim cWz= Ci , (i= 1, 2, n), 

ft =00 

ana* if (yi, ■ ■ • y n ) an ^> {yVy • ■ • • Vn?) are solutions of (l) ana" (16) respectively 
which satisfy the initial conditions : 

yi (c) = d, V? ] (o) = c?\ (i = 1, 2 n) , 

^e» tfAe functions yf ) {x) approach uniformly the functions yi(x) as k becomes 
infinite. 

* It should be noticed that we need not here exclude the possibility that the number of discontinu- 
ities of some or all of the functions a$>, $*> become infinite with ft. Cf. the special case mentioned 
below. 
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If, in particular, the coefficients of (1) are continuous throughout (I), we 
may form the equations (16) by dividing the interval (I) into tn k portions, the 
integer ni k being supposed to increase without limit as h becomes infinite, and all 
portions into which (/) is divided being supposed to approach zero as their limits. 
Let us further denote by xf\ xf\ . . . . x^ % points lying one in each of these sub- 
intervals ; and let us define af and $*> throughout the Z th of these intervals by 
the formulae : 

af{x) = a i} {xf), f(*)=AW. 

The coefficients of (16) thus consist, for each value of h, of a succession of 
constant values, so that the functions yf> can be expressed in terms of elementary 
functions. The method here sketched for approximating to the solution of (1) 
was developed for a special case at considerable length by Radakovic.* 

The results of the present paper can be applied without difficulty to the case 
of a single linear differential equation of the n th order, since such an equation is 
equivalent, as is well known, to a system of the form (l).f 

Harvard University, Cambridge, Mass. 

* Monatshefte flir Math. u. Physik, vol. 5 (1894), p. 193. 

f The results thus obtained -will include, among others, the theorems established by the writer in 
the Bulletin of the American Mathematical Society for March, 1899, on p. 276. 



